Proton Single Particle Energy Shifts due to Coulomb Correlations 
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A theoretically consistent approach to the calculation of the Coulomb cor- 
Q^". relation corrections to the single-particle energies is presented. New contribu- 

tions to the single-particle energies previously overlooked in the literature are 
now identified and taken into account. We show that the interplay between 
. the Coulomb interaction and the strong interaction, which is enhanced in the 

<n : 

nuclear surface, leads to an upward shift of the proton single-particle levels. 
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This shift affects the position of the calculated proton drip line, a shift towards 
decreasing Z. We describe briefly a similar mechanism which is at work for 



neutron levels. The same mechanism is responsible for significant corrections 
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to the mass difference of the mirror nuclei (Nolen-Schiffer anomaly) and to 
the effective proton mass. 



PACS: 21.10.Sf Coulomb energies — 21.10.Dr Binding energies — 21.10.-k Nuclear energy 
levels 

The main part of the Coulomb energy in nuclei is given by the Hartree contribution 
and to a reasonable accuracy this can be computed as the energy of a uniformly charged 
sphere and is thus proportional to Z 2 /A 1 ^ 3 . There are a number of corrections, some of 
them rather subtle, arising from the interplay between the Coulomb and nuclear forces. The 
Nolen-Schiffer anomaly in the binding energy of mirror nuclei ffl is a case in point. In Ref. 
H we have shown that a specific many-body mechanism (outlined concisely below) leads to 
an enhancement of the Coulomb energy in the nuclear surface region and in this way one can 
account for a major part of this anomaly. This effect results in a systematic contribution to 
the nuclear binding energy, which scales as oc Z 2 ! 3 . We are going to demonstrate that the 
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mechanism should be taken into account when calculating a number of nuclear properties. 
In this Letter we study the role of this new many-body effect on the single-particle proton 
energy levels, on the location of the proton drip line and on the proton effective mass. 

We shall operate within the density functional theory The ground state energy of 

a nucleus E is given by a sum of two functionals (in the absence of pairing correlations): 



where the symmetric part Fo[p n (r), p u (r)] = Fo\p u (r) , p n (r)] is due to the (strong) isospin 
conserving nuclear forces, while Fc ou i is due to the (weak) Coulomb interaction. We shall 
neglect in our analysis several easy to include terms: the trivial contribution in the kinetic 
energy, arising from proton-neutron mass difference, the contribution of the Charge Sym- 
metry Breaking (CSB) forces flj]]. For the sake of simplicity of the presentation, we shall not 
display explicitly the spin degrees of freedom and the contribution arising from the spin- 
orbit interaction, even though we have included them in the actual calculations. The proton 
and neutron densities are defined as usual 



where n^ii^i and n v i^(j) v i are proton and neutron quasiparticle occupation numbers and 
single-particle wave functions respectively. The well known Skyrme functional, see e.g. Refs. 
|||9|], can be considered as one possible realization of F Q . The Coulomb energy functional is 
given by: 
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F [p n (r),p u (v)] + F C oui[pn(r),p u (r)}, 
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with p 7r7r (ri, r 2 ) being the exact two-proton density distribution function. Xwn( r i, r 2, ioS) and 
X°( r i> r 2 ; i^) are the full and free proton linear response functions respectively, evaluated at 
the imaginary frequency iu. Eqs. ([5]) and (|6|) represent the Hartree and Fock contributions 
to the Coulomb energy. The exchange term is written here in a somewhat unusual way [0, 
through the linear response function of the noninteracting protons, %° , since upon integrating 
along the real axis of the complex uj plane one has 

/ — Imx° Oi, r 2 , lo) = [p^ri, r 2 )p n (r 2 , r x ) - <J(n - r 2 )p 7r (ri)], (8) 
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where /^(ri,^) is the proton single-particle density matrix |10 |. In Eqs. (|6],[7]) we have 



performed a Wick rotation in order to evaluate the integral along the imaginary axis. Upon 
taking into account the (strong) residual interaction, the linear free response function xl 
should be replaced with the full response function Xnn and one thus readily obtains the 
expression for the Coulomb correlation energy (|7|). By considering all three contributions 
we therefore account for all diagrams in first order in e 2 , in the (weak) Coulomb 
interaction. The (strong) nuclear interaction is treated to all orders. We are going to 
concentrate on a study of the contribution to the single particle energies and the effective 
mass due to Eq. (|7|). In passing we remark that there are no ambiguities with double 



counting, rearrangement energy, core-polarization correction and so forth, see Refs. [|11[] for 
a discussion of this issues, within the density functional formalism. 

By expressing through the linear response function we can easily calculate a num- 
ber of functional derivatives needed below. At the same time one can clearly disentangle the 
contribution of various modes to this part of the energy density functional as well [|J . In 
Ref. we have shown that the main contribution to comes from the surface collective 
isoscalar excitations. In the isoscalar channel the particle-hole residual interaction has a 
strong density dependence, changing from a relatively weak one inside nuclei to a strong 



attractive one in the surface region |TS| . Because of the attractive character of the isoscalar 
residual particle-hole interaction, the nuclear surface (where the matter density is low) is 
very close to instability |§. Low density homogeneous nuclear matter is manifestly unstable 



and this shows itself in the fact that the linear response function has a pole for imaginary 
values of to ||. In finite nuclei and semi-infinite nuclear matter this singularity is smoothed 
out and becomes a prominent peak of the response function at low frequencies in the surface 
region. This proximity of the nuclear surface to instability is the reason why the contribu- 
tion of the (weak) Coulomb exchange interaction is strongly enhanced and the correlation 
Coulomb "correction" to the mass formula has as a result a predominantly surface character, 
i.e. Fg% oc Z 2 /\ 

Let us turn now to the calculation of the proton single-particle energy corrections due 



to the presence of Fq£L in the energy functional. Using Landau's variational equation [16 



£*i=x , 9 

one obtains for the proton single-particle energy shift Ae^ i the following expression 

e 2 r dv 1 dr 2 du 8[x^{vu «w) - x°(r i; r 2 , iu)] 
2 J 27r|r! - r 2 | dn^i 

The variational derivative 5x^ / Sn n i has the simple functional form, 

6£(r lt T 2 ,uj) = ^ + £jt } + (w _ ( ri )^ Ao(r2 ), (ii) 

with G n (ri, r 2 , u) being the single-particle proton propagator in the selfconsistent nuclear 
potential. The linear response function Xmr is obtained by solving the usual matrix functional 
equation (Landau zero sound or RPA) 

Xij = X%j + X°iRikXkj, (12) 

k=TT,l/ 

where Rik is the irreducible (strong) particle-hole interaction and i,j and k stand for the 
isospin variables, tt and v for protons and neutrons respectively. From this equation one 
derives the following equation for 5Xi™ / fin* i 



^Xij _ $Xi j 



S Xm p , ^p SRjk On S Xkj 

tiikXkj + Xi 7 Xkj + Xi^ik- 



(13) 



The operator solution can be easily obtained and one thus can show that the relevant 
quantity in the integrand of Eq. (0) has the following structure 
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where for the sake of clarity we have suppressed the isospin, spin and spatial coordinates 
and the corresponding summations and integrations. 







FIG. 1. The diagrammatic representation of corrections to the single-particle nucleon propaga- 
tor, corresponding to the Coulomb correlation energy discussed in this Letter. The solid oriented 
lines represent the nucleon propagators, the dashed lines the Coulomb interaction, the thin wavy 
line the strong nuclear interaction it^ri, r2) and the double wavy lines stand for the "screened" 
strong nuclear interaction, obtained by solving the diagrammatic equation shown in the third row. 

In Fig. 1 we show the diagrammatic representation of the corrections to the single- 
particle nucleon propagator, corresponding to the addition of the Coulomb correlation energy 
to the nuclear density functional given by Eq. (7). Strictly speaking, we have omitted in this 
diagrammatic representation the contribution arising from the last term in Eq. flip), which 
depends on the functional derivative of the irreducible particle-hole interaction 5R/5n n i. 
We shall discuss in more detail the relevance of such a correction in a future publication. 
To the best of our knowledge this correction was not explicitly discussed in literature in 
this context. A superficial analysis would suggest that this type of correction is perhaps 
small. The two diagrams shown in the upper row were considered before by other authors 
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ljj , however the third diagram, shown in the middle row, was, surprisingly, overlooked. 
All three diagrams are of the same order in the electromagnetic coupling constant e 2 . Only 
when considering all three diagrams does the correction to the single-particle energy satisfy 
the natural Landau's variational equation (9). Besides being theoretically consistent, the 
present approach has a new physical feature as well: there are corrections to the single- 
particle neutron properties too. While for protons all three diagrams have to be considered, 
for neutrons there is a Coulomb correlation correction to the single-particle energies arising 
from the third diagram alone. 



We do not consider in this Letter the so called AKW diagrams [12|, which are automat 



ically taken into account in selfconsistent calculations of nuclear masses @. Auerbach [13 



has shown that the corrections corresponding to these diagrams are too small to account for 
the Nolen-Shiffer anomaly. 

In Ref. we have used a simple separable model for the residual interaction Rik 

Rik{*i, r 2 ) = A — — — — d(S2i - S2 2 ), (15) 

where Vi(r) is the proton/neutron single-particle potential in a spherical nucleus and S(Qi — 
Q2) is a Dhac function in angle variables. A is chosen so that the dipole linear response has a 
pole at uj = 0, corresponding to the spurious mode. This type of residual interaction has been 
studied extensively and it leads to a satisfactory description of nuclear collective 

modes. We have used the same nonselfconsistent approach described in detail in Ref. [Q] in 
order to estimate the magnitude of the correction Ae^i- For various single-particle proton 
levels around the Fermi level and the proton threshold, the calculated shifts Ae n \ are in the 
interval 0.1 — 0.3 MeV in light (A = 16) and medium (A = 40 — 48) nuclei. In performing 
these calculations in each nucleus we have included collective modes with multipolarities 
up to ss 2A 1 / 3 and energies up to « 225 MeV, see Ref. 0. Both in multipolarities and in 
energy convergence was attained in our calculations. The values of Ae w 1 thus obtained are 
of the same magnitude as the Nolen-Schiffer anomaly. 

It is instructive to cross check these results using an independent approach. A new type 
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of nuclear density functionals has been recently introduced in Refs. PH, see also Ref. ||. 
The main reason for seeking new functionals is to obtain a significantly more accurate repro- 
duction of the nuclear properties (bindings energies and matter distribution of finite nuclei 
and infinite neutron and symmetric nuclear matter properties over a wide range of densities 
simultaneously) than one can achieve with the plethora of existing density functionals. A 
key ingredient was the introduction of the Coulomb correlation energy contribution, along 
the lines suggested earlier by us in Refs. 0. However, while we have presented arguments 
for a significant surface contribution into the nuclear Coulomb correlation energy, Fayans 
H has chosen to parametrize the Coulomb correlation energy as a volume term: 

(16) 

Here po — 0.16 fm -3 , a = 1/3 and from a fit of the masses and radii of 100 medium and heavy 
spherical nuclei Fayans determined hc ou i = 0.941, see also Ref. 0. The Coulomb correlation 
energy determined by Fayans thus almost exactly cancels the Coulomb exchange energy, 
which in the Slater approximation is formally given by the same formula with hc ou i — — 1 
and (7 = 0, namely 

A quick estimate of the Coulomb correlation energy contribution shows that 

FcZi « \ (^) V3 e 2 h Co ui J drp w (r) » 0.4 ZMeV, (18) 

which leads to a typical shift Ae n \ ~ 0.4MeV of the same order of magnitude as estimated 
by us independently. An upward shift of this magnitude of the last occupied proton level 
in a nucleus near the proton drip line is equivalent to a shift of the proton drip line in the 
direction of decreasing Z by a few units. A shift of this magnitude for a neutron level would 
be equivalent to changing the mass number by up to 5 units, see for example |19[ . A similar 
shift due to the Coulomb interaction arises for neutron levels, but we shall not discuss this 
mechanism here (see however Fig. 1 and the ensuing discussion). 
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In a recent preprint Brown et al. ]14| show that one can generate essentially an innumer- 



able range of phenomenological nuclear density functionals of the Skyrme type, with which 
the Nolen-Shiffer anomaly can be accounted for, while at the same time nuclear masses and 
radii can be calculated selfconsistently with very good accuracy. This lack of uniqueness 
for the nuclear density functional is due to the lack of an uncontroversial theoretical under- 
pinning. It is hoped that the present results will bring the goal of a complete theoretical 
understanding within our grasp. 

One can show that there is another related effect due to the Coulomb correlation energy, 
a noticeable renormalization of the proton effective mass near the Fermi surface. In the case 
of homogeneous nuclear matter the effective mass is defined as [plj 

1 1 de(p) 



(19) 

P=PF 



m* pp dp 

where pf is the Fermi momentum. One can show that the effective mass renormalization 



Am can be obtained from the relation 

Am e 2 d [ f dqdu Sxoiq,^) 1 



(20) 



m*(m* + Am) pp dp [J (2tt) 3 5n{p) q 2 [l — R(q, iu)xo(q, iui)] 2 J ' 

where m* is the proton effective mass computed in the absence of the Coulomb correlation 
energy and n(p) is the quasiparticle occupation number probability of the state with linear 
momentum p. Using the following approximate formula (which becomes an identity at the 
point where the compressibility is vanishing) 



/'-;<;- P p 



-— KVf- |p + q|)5Mp- (p + q), (21) 



dp 6n(p) 

one can reduce the above expression for the effective mass shift to 



1 1 e 2 r l xdx , nn . 



m* + Am m* 2np F J-i [1 — R(q(x), 0)xo(q( x )i 



where q(x) = pp\J 2(1 — x). At the point where the compressibility tends to zero, the de- 
nominator [1 — R(q(x), 0)xo(q(x), 0)] vanishes at x = 1 (q — 0). We have shown in Ref. ]2[] 
that the nuclear surface of finite nuclei is rather close to this regime. Since the integrand in 
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Eq. fl22|) is positive, the integral diverges and the effective mass thus vanishes in an infinite 
homogeneous system. In finite nuclei this divergence is smoothed out and the pole singu- 
larity becomes a narrow surface peak ||. The net effect is that the proton effective mass 
becomes smaller than the neutron effective mass, and therefore the proton level density is on 
average smaller than the neutron level density in nuclei. B.A. Brown argues that a somewhat 
similar in medium effective nucleon mass renormalization can account for the Nolen-Schiffer 
anomaly ||. An earlier QCD sum-rule approach [^] substantiate such a claim. This effect, 
as well as the contributions arising from CSB forces @, lead mainly to volume effects, while 
the many-body mechanism discussed by us leads to a surface contribution. 

There is thus hope not only to generate in the near future extremely accurate nuclear 
density functionals, but also to be able to understand the nature of various rather subtle 
contributions. We anticipate as well that the study of the Nolen-Schiffer anomaly within 
the framework of the density functional approach will allow us to estimate the magnitude 
and the character (volume versus surface) of the CSB many-body effects. 
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